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Abstract 

The Oklo data constrain the depth of the nuclear potential well to a 
small margin of error, determined by various unknowns, such as the reac- 
tor temperature. However, we show that even these small variations could 
allow the U-238 half-life to vary by more than one order of magnitude. 
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1 Introduction 



Beginning with Shlyakhter [1] in 1976, Oklo natural reactor data have been 
used by many authors [2-6] to put limits on the time variation of " constants" 
such as the fine structure constant and the strong coupling constant. As noted 
by Damour and Dyson [2] , such constraints generally refer to the value of the 
"constants" at the time of the reactor shut down, since most prior evidence is 
diminished by the reaction process. Also, it cannot be ruled out that two or 
more factors may change concurrently, and one may cancel the effect of the 
other on whatever we are able to calculate for comparison with experiment. 

In this work, we will use the constraints provided by the samarium- 149 cross 
section to provide a self-consistent limit on the value of the uranium-238 decay 
constant at the time of reactor shutdown. This implies a limit on possible differ- 
ences between the strong coupling constant at the time of reactor shutdown and 
today. Pierronne and Marquez [7,8] developed the complex eigenvalue method 
for calculating alpha-decay half-lives by matching the wavefunction of the alpha 
particle, and its derivative, across a boundary between a square well represent- 
ing the nuclear interior and a coulomb potential representing the exterior of 
the nucleus outside the range of nuclear forces. Chaffin, Banks, Gothard, and 
Tuttle [9,10] have further developed this method into Mathematica notebooks 
permitting easy calculation of relevant variations. 

Experimental data on alpha particle scattering [11-14] have long been known 
to validate a cluster model in which an alpha particle moves in the average 
field provided by the rest of the nucleus. Four-nucleon transfer reactions were 
observed on target nuclei that are the daughters of an alpha emitter. The 
comparison of reaction rates with the corresponding alpha-decay rates indicates 
the tendency for the four transferred particles to form a true alpha particle. 
The data give evidence that nodes in the alpha particle wavefunction do in fact 
exist. When the nuclear potential seen by the alpha particle is represented by a 
simple square well, the effect of changing the depth of the well is to change the 
alpha particle wavefunction. Usually this change is slight, but at certain critical 
values the number of nodes in the wavefunction can change precipitously, with a 
corresponding change of more than an order of magnitude in the decay constant. 
Hence, we find that, although the Oklo data provide Shlyakhter's very stringent 
limits on the variation in, for instance, the Sm-149 absorption cross-section for 
neutrons, this may not translate into such stringent limits on nuclear alpha- 
decay half-lives. 

In section 2 of this paper we will present the method of calculation of the 
Samarium cross-section at the time of the Oklo reactions, and indicate how 
this may be used to find the depth of the square well representing the nuclear 
potential. Then in section 3 we review the Pierronne and Marquez [8] method for 
calculating the decay constant from the depth of the well and other parameters. 
We then present the results of a calculation for the case of alpha-decay of U-238. 
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2 Using Oklo Data to find the Samarium-149 
Cross-Section 



We have used actual data of element concentrations from Oklo samples to find 
the implied Sm-149 absorption cross-section for neutrons. Procedures very sim- 
ilar to those of Damour and Dyson [2] and Fujii et al.. [3] were used. As Fujii et 
al. explained in their Appendix B, coupled linear differential equations may be 
written down for the time derivatives of the Nd-143, Sm-147, U-235, and U-238 
concentrations. These equations may be solved for the restitution factor C rz , 
epi-thermal index r rz , and the neutron fluence appropriate for each sample. 
For each sample the ratio k = N® 35 / 'N nat , the ratio of the initial U-235 abun- 
dance to the natural abundance of samarium, is also found using Fujii et al. 's 
procedure. 

As explained in the Damour and Dyson [2] and/or Fujii et al. [3] papers, 
coupled linear differential equations may also be written for the Sm-147, Sm- 
148, and Sm-149 concentrations. To obtain a working equation we use equation 
(16) of Fujii et al. [3], corrected for a typographical error in the subscript for 
the Sm-149 absorption cross-section, yielding the result for the abundance ratios 
at the end of the reactor activity. If the time for the end of reactor activity is 
ti, and N147 , TV^s, and ./V149 are the concentrations of Sm-147, Sm-148, and 
Sm-149, respectively, the appropriate equation is: 



Here the R's are the relative fractional natural abundances of the isotopes of 
samarium indicated by the subscripts, Y147 and Y149 are the fission yields of 
Sm-147 and Sm-149, respectively. Using this equation, we used an iteration 
procedure called FindRoot in Mathematica to obtain the samarium-149 absorp- 
tion cross-section. 

Results are given in Table I, where they are compared to Fujii et al. 's 
results. 

The table gives the cross-sections obtained for samples inside the core of 
reactor zone 10, and average cross-section of 90.29 ± 7.41 kilobarns. The errors 
are the one standard deviation statistical errors. The temperature of the Oklo 
ores at the time of the reactions is unknown but thought to have been around 
300 to 400 degrees Celsius. This temperature would, of course have varied over 
the lifetime of the nuclear reactor, and is an important systematic source of 
error in determining, from the data, the Sm-149 cross-section at the time of the 
reactions. 



ATi4 9 (ti) 




(1) 



CT /235 y 149 
CTa-o- 1 49 



[exp(-CTi4g0ti)-cxp(-CT a 0ii)]+R™49 exp(-(Ti49^ti) 
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3 Finding the U-238 Decay Constant 



In a simple nuclear model in which a square well is used for the potential seen 
by the alpha particle inside the nucleus, the Sm-149 cross-section is related to 
the depth of the potential well. The cross-section as a function of energy is 
given by the Brcit-Wigncr shape: 

r r 

" 7 (2) 



{E-E r r+Tl t /A 

In this expression the widths T do not depend directly on the depth Vo of the 
potential well, but the resonant energy E r is a constant minus Vo , as was shown, 
for example, by Weber, Hammer, and Zidcll [15]. The spread in values of the 
samarium cross-section found in section 2, thus can be related to a spread in 
values of Vo- For U-238, the half-life of 4.47 x 10 9 years is reproduced by a Vo 
of 110.5 Mev. From the data of section 2, the integrated Sm-149 cross section 
is found to be 90.3 ± 7.4 kilobarns (see Table I). The samarium resonant energy 
is E r =0.0958325 eV and the corresponding one standard deviation spread is 
±0.0102 eV. While this spread might seem small, the small change in Vb can 
lead to a large change in the alpha-decay half life. 

To test the variability of the decay constant, Mathematica notebooks were 
written. For the square well potential on the inside of the nucleus, the Mathe- 
matica notebook gave answers essentially equivalent to those of Chaffin, Gothard 
and Tuttle [10], which used a harmonic oscillator potential for the interior re- 
gion, where the nuclear potential is felt by the alpha particle. In the course of 
this work, it was discovered that, as the nuclear potential well depth is changed, 
and the nuclear radius changes slightly, it is possible to have a sudden change in 
the number of nodes of the real part of the alpha particle's wavefunction. This 
was modeled for both the harmonic oscillator and square well potentials, with 
nearly the same results for either notebook. 

The change in the number of nodes causes the probability of tunneling to 
change by about a factor of 13.5, as shown by the discontinuity in the graph of 
Figure 2. The graph shows a large, precipitous rise in the value of the decay 
constant A, at certain critical values of Vo, similar to the postage-stamp function. 



4 Conclusions 

The calculations of the decay constant versus well depth Vo show that even a 
small changes in the strong coupling constant could result in a U-238 decay con- 
stant change by a factor of 13.5, more than one order of magnitude. The effect 
is tied to the change in the number of nodes in the alpha-particle wavefunction 
when the well depth passes critical values. 
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Figure Captions 

Figure 1. Sudden change in the number of nodes. The harmonic oscillator 
wave function for well depths of 58 MeV (a) and 54 MeV (b) . The x-axis is the 
radial coordinate of the alpha particle, T — pj (2ry), where p and rj arc defined in 
Green and Lee [16]. Figure la shows the harmonic oscillator wave function for 
a well depth of 58 MeV. Figure lb shows what happens when the well depth is 
changed to 54 MeV, without changing the alpha particle energy. If one counts 
the number of nodes in Figure la, there are nine, not counting the ones at zero 
and infinity. For Figure lb, there are eight nodes, a reduction by one. 



Figure 2: The Decay Constant A of U-238 versus Well Depth. 
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Sample 


(Ti49 (kb) 


Fujii et al. results 


SF84-1492 


98.1 


99.0 


SF84-1485 


83.0 


83.8 


SF84-1480 


95.0 


96.5 


SF84-1469 


85.0 


85.6 


standard deviation 


7.41 


7.64 


average 


90.3 


91.2 



Table 1: Sm-149 Effective Cross Sections 
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Figure 1. Sudden change in the number of nodes. The harmonic oscillator 
wave function for well depths of 58 MeV (a) and 54 MeV (b). The x-axis is 
the radial coordinate of the alpha particle, T= p/(2r|), where p and r| are 
defined in Green and Lee [16]. Figure la shows the harmonic oscillator 
wave function for a well depth of 58 MeV. Figure lb shows what happens 
when the well depth is changed to 54 MeV, without changing the alpha 
particle energy. If one counts the number of nodes in Figure la, there are 
nine, not counting the ones at zero and infinity. For Figure lb, there are eight 
nodes, a reduction by one. 
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Figure 2: The Decay Constant X of U-238 versus Well 
Depth . 



